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We study inelastic energy relaxation in graphene for low energies to find out how electrons scatter 
with acoustic phonons and other electrons. By coupling the graphene to superconductors, we create a 
strong dependence of the measured signal, i.e., critical Josephson current, on the electron population 
on different energy states. Since the relative population of high- and low-energy states is determined 
by the inelastic scattering processes, the critical current becomes an effective probe for their strength. 
We argue that the electron-electron interaction is the dominant relaxation method and, in our 
model of two-dimensional electron-electron scattering, we find a scattering time r e _ e = 5 . . . 13 ps 
at T = 500 mK, 1-2 orders of magnitude smaller than predicted by theory. 

PACS numbers: 72.80.Vp, 74.45.+C, 72.15.Lh 



I. INTRODUCTION 

The special electronic structure of graphene shows up 
in its electronic properties ^ Most attention has been 
paid to the electronic conductivity^ which, due to the 
strong energy dependence of the density of states, can 
be tuned significantly with a gate voltage. For dirty 
graphene, this gate dependence is furthermore modified 
by elastic scattering^ due to potential inhomogenities 
forming charge puddles. At high voltages, also inelas- 
tic scattering due to optical phonons appears^, but the 
low-energy inelastic scattering due to electron-electron 
(e-e)&£ or electron- acoustic phonon scattering (e-ph)^ 9 - 
does not directly influence the conductivity because the 
mean free paths for them are typically larger than the 
elastic mean free pathi^ 

In this paper we study the effect of low-energy inelastic 
scattering in graphene by using Josephson critical current 
measurements to determine heat transport in the system. 
The idea is to apply a heater voltage and to measure the 
increased temperature via nearby thermometers. The 
measurement is performed at sub-kelvin temperatures 
and low voltages, thereby providing access to low-energy 
inelastic scattering processes and allowing to disregard 
scattering from optical phonons. To perform the mea- 
surement we use three superconducting electrodes fab- 
ricated on graphene. Two of them (thermometer elec- 
trodes, C and R in Fig. [TJ lie close to each other, so that 
we can measure a finite supercurrent through them, and 
the third one (heater electrode, L in Fig. [T]) is used for 
heating the system. The supercurrent is sensitive to the 
electron temperature or, more accurately, to the electron 
distribution function on the graphene region between the 
superconductors^ and therefore acts as an electron ther- 
mometer. 

We make two measurements: First we measure the 
supercurrent as a function of temperature, and thus cal- 
ibrate the thermometer and find parameters for a micro- 
scopic theory describing the supercurrent in the junction. 
Second, we apply a voltage to the heater electrode, sup- 
plying heat into the electron system, and measure again 
the supercurrent in the presence of the heater voltage. 
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FIG. 1. (color online), (a) Optical image of the studied 
graphene flake with three Ti/Al electrodes marked L, C and 
R. The contour of the flake is highlighted by a a dashed line, 
(b) Differential resistance of the long (L-C) and short (C-R) 
sections as a function of gate voltage at T = 4 K. 



The magnitude of the supercurrent in the presence of 
the heater voltage is sensitive to the strength of inelastic 
relaxation inside graphene, allowing us to measure it. 

The Joule heat generated in the presence of a 
bias voltage is dissipated either to the electrodes or 
to phononsi 8 ' 9 ' 20 However, since the electron-acoustic 
phonon coupling is so weak in graphene, most of the heat 
escapes into the electrodes, even though this process is 
blocked at low energies by the superconducting gap A. 
The escape is possible as a result of processes transferring 
excitations from low energies to above the gap, in par- 
ticular multiple Andreev reflections^ and e-e scattering. 
The former tends to broaden the electron distribution be- 
low the gap and thus to increase the effective temperature 
at these energies, while the latter drives the distribution 
towards a quasiequilibrium (Fermi function) form hav- 
ing a lower effective temperature. Our thermometer is 
most sensitive to the distribution at energies below A, 
and therefore it is a sensitive probe of the e-e scattering 
strength in graphene. But to extract the magnitude of 
the e-e scattering we have to abandon the simple effective 
temperature description used for example in Ref. [ID and 
rather solve a full kinetic Boltzmann equation, with e-e 
scattering included explicitly with a collision integral. 

This paper is organized as follows: In Sec. [Til wc 
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describe the experimental setup used to carry out the 
measurements. In Sec. IIIII we formulate our theoreti- 
cal model and consider the different sources of inelastic 
relaxation. In Sec. II VI we combine the theoretical and ex- 
perimental results to provide an estimate for the strength 
of relaxation in the system. Finally, we discuss the im- 
plications of these results in Sec. [V] 



II. EXPERIMENTS 

An optical image of the studied monolayer graphene 
sample is shown in Fig. [Ha). The 2.8 /itm long and 4.0 
/xm wide graphene area in between leads L and R is par- 
tially interrupted by a 1.0 fjm wide lead C. This latter 
lead is 1.7 ^m and 0.4 ^m far from leads L and R, re- 
spectively, with all the distances measured between leads' 
internal edges. The graphene flake has been exfoliated 
with a semiconductor wafer dicing tape and deposited 
on top of a 250 nm thick SiC>2 layer. The oxide isolates 
the graphene flake from a highly p-doped Si substrate 
used as a back gate in our experiments. Three Ti/Al 
(10 nm/50 nm) metallic contacts were patterned by us- 
ing standard electron-beam lithography techniques, and 
evaporated in ultra-high vacuum. A 10 -10 mbar vacuum 
during the metal evaporation guarantees highly trans- 
parent contacts, which are needed to observe proximity- 
induced supercurrents. The sample was measured at low 
temperatures, down to 80 mK, in a dry dilution cryostat 
BF-SD250 from Blucfors. The sample contacts are elec- 
trically connected to room temperature electronics via 
one-meter long thermocoaxes, low-pass RC filters (cut- 
off frequency of 1 kHz) and one-meter twisted pairs, pro- 
tecting the sample from the room-temperature electri- 
cal noise. Below the critical temperature of ~ 600 mK, 
the Ti/Al leads become superconducting, resulting in the 
formation of three different superconductor-graphene- 
supcrconductor (SGS) junctions. 

Figure QJb) shows the gate voltage (V g ) dependence 
of the differential resistance (Rd) of the graphene sec- 
tions L-C, C-R and of the whole flake (L-R), at 4 K. The 
entire section L-R presents a peak in Rd at a gate volt- 
age Vcnp = — 11 V. This resistance peak is associated 
with a minimum charge carrier density and takes place 
at the charge neutrality point (CNP). The negative value 
of Vcnp indicates that the flake is n-doped in the absence 
of gate voltage. The change in the resistance is smaller in 
the p-doped region than in the n-doped one. This asym- 
metry comes from the n-doping by the Ti/Al leads^~— 
and the resulting formation of p-n junctions when bulk 
of the graphene is p-doped by the gate for V g < Vcnp- 
For V g > Vcnp , we estimate from the Drude model a mo- 
bility of 3500 cm 2 V s -1 and a mean-free-path of 70 nm 
at V g = 30 V. The electrical transport through the stud- 
ied sections is therefore diffusive. The resistance Rn of 
the short section (C-R) is continuously decreasing with 
V g and no resistance peak is observed within the investi- 
gated gate voltage range. The short section is thus always 
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FIG. 2. (color online). Current-voltage characteristic of the 
short section C-R at V g = —30 V obtained by successively 
sweeping up (blue curve) and down (red curve) the bias cur- 
rent at 80 mK. The SGS junction jumps from the zero- voltage 
state into the resistive state at the critical current I c and 
comes back into its initial state at the retrapping current I r . 
Inset: IV i 
T=80 mK. 



n-doped. This is because the contacts affect the sample 
on a scale of micrometers so that the average doping in 
the short section is stronger. 

We now focus on the current-voltage (/V) characteris- 
tics of the long and short sections measured at 80 mK for 
different gate voltages. The IV curves generally contain a 
supercurrent branch characterized by a zero- voltage state 
when the bias current is kept below the critical current 
I c (see Fig. [2]). At J c , the SGS junction jumps into the 
resistive state. The initial superconducting state is re- 
covered when the current is biased below the so-called 
retrapping current I r . For the short sample, at 80 mK, 
I r always differ from I c leading to a hystcrctic IV curve. 
The critical current I c depends on the gate voltage and 
decreases when the normal resistance increases as seen 
in Fig. [T] When the gate voltage is tuned from —30 to 
30 V, the critical current increases from 47 to 120 nA 
and the normal resistance is lowered from 708 to 337 fi, 
respectively. 

Fig. [3] shows the current-voltage characteristics of the 
sections L-C and L-R at V g — —3.4 V. Both IV curves 
present a supercurrent branch with the same critical cur- 
rent of 16 nA, and identical differential resistance values 
(650 f2) at sufficiently low bias. This is understood by 
the presence of a supercurrent through the short section 
(C-R) keeping leads C and R at the same potential. The 
transition of the short sample into the resistive state is 
identified by a second voltage jump. This jump takes 
place at a critical current of 156 nA for the section L-C 
and 104 nA for the section L-R. The differential resis- 
tances increase to 685 and 810 SI for sections L-C and 
L-R, respectively. 

The IV characteristics of our SGS junctions are 
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FIG. 3. (color online). IV curves of sections L-R and L-C at 
V g = —3.4 V and T = 80 mK. Three voltage transitions are 
highlighted by the arrows. 




FIG. 5. (color online). Critical (circle points) and retrap- 
ping (cross points) currents as a function of the voltage VL-c 
across the section L-C. The triangle points corresponds to the 
deduced electronic temperature. 
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FIG. 4. (color online). Temperature dependence of the 
critical (solid line) and retrapping (dashed line) currents at 
V g = —3.4 V. Inset: Current-voltage characteristic of the 
short section C-R at V g = —3.4 V for different temperatures. 
The current is swept up from -120 to 120 nA. 



Keeping the bath temperature at 80 mK, the electronic 
temperature can be changed by injecting a dissipativc 
current in between leads L and C. As shown in Fig. [5] 
for V g = —3.4 V, the critical current decreases in the 
short section (C-R) as the voltage VL-C across the long 
section L-C goes up. The retrapping current remains ini- 
tially almost constant at low voltages and decreases at 
around 0.15 mV. By assuming that the electronic dis- 
tribution follows a Fermi distribution, we can directly 
relate the critical current to the electonic temperature. 
We find that the temperature amounts to around 0.52 K 
at Vl-c = 0.35 mV. However, as we show in Scc. lIII| the 
electronic distribution function may differ from a Fermi 
distribution in which case the electronic temperature is 
not properly defined. Consequently, the electronic tem- 
perature directly deduced from the measurement of I c is 
only an effective temperature (see Eq. below). 



III. THEORETICAL MODEL 



strongly changing with temperature. Fig. [4] shows the 
temperature dependence of the TV curve of the short sec- 
tion (C-R) at V g = —3.4 V. The critical current decreases 
with increasing the temperature. The retrapping current 
remains almost constant until 0.3 K and then goes down. 
The hysteresis of the IV curve is reduced as the temper- 
ature goes up and disappears at 0.45 K. Similar results 
are found at different gate voltages, and also in the long 
section L-C. The temperature dependence of the critical 
current is used in Sec. IIVI to extract the charge carrier 
mean free path. From the shape of this dependence we 
can already tell that we are in the long junction regime, 
where the length of the junction L is longer than the su- 
perconducting coherence length £. Most importantly, the 
strong temperature dependence of the critical current al- 
lows us to use SGS junctions as electronic thermometers. 



We compose a model where a graphene flake is di- 
vided into two parts by superconducting electrodes with 
energy gaps A and treat the system as effectively one- 
dimensional with the essential dimension aligned along 
the x-direction. This setup is illustrated in Fig. |6] where 
the two distinct regions are numbered as 1 and 2. 



A. Distribution function and supercurrent 

Physical observables, such as the supercurrent, can 
now be determined from the electron distribution func- 
tion /(e) which is a function of both energy e and posi- 
tion x, although below the latter is not explicitly writ- 
ten down. The distribution function satisfies the time- 
independent diffusion equation (disregarding the proxim- 
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ity effect) 

- D ^ =I{f{€)) - (1) 

Here D is the diffusion constant, related to the Fermi 
velocity vp = 10 6 m/s and the transport relaxation time 
r for elastic scattering by 

D = v%t/2, (2) 

and /(/) is the collision integral for the inelastic processes 
in the flake (most importantly e-e and e-ph scattering). 
At the graphene-supcrconductor interfaces we have two 




FIG. 6. (color online). Sketch of the model with graphene 
flake divided into two regions: the heater (region 1) and ther- 
mometer (region 2). The relative sizes of the regions are de- 
termined from critical-current measurements as explained in 
the text. 

kinds of boundary conditions depending on whether we 
are at the end points x — and x = L or at the boundary 
of regions 1 and 2, at x — x c . At the end points^ 

m = f L/R (e), \e\>A (3a) 
/( At + e ) = l-/( At - e ), | e |<A (3b) 

d x f(ti+e)=d x f(»-e), |e|<A. (3c) 

Here, /i/ij(e) is the equilibrium Fermi function at the 
corresponding potential, i.e., /l(c) = fo(e — eV) and 
fn(e) = /o(e) as / (e) = [cxp [(e - M )/fc s T bath ] + If 1 
at base temperature Tb a th and chemical potential fj, with 
respect to the Dirac point. Below the gap, these bound- 
ary conditions conserve the balance of positive and neg- 
ative charge excitations (|3b[) and ensure that there is no 
energy current entering the superconductor (|3c[) . due to 
Andreev reflection. Above the gap, we simply have conti- 
nuity of the particle distribution across the contact since 
the interface is assumed transparent (|3a[) . At the divid- 
ing superconductor at x — x c , we require 

/(*)!*=*._ = mi= Xa+ (4a) 



/(e) = / (e), |e|>A (4b) 
/( M + e ) = l-/( M -e), |e|<A (4c) 

d* ffa + e)\ x=Xe+ - d x /( M - e)\ x=Xo+ , \e\ < A. 

(4d) 

In addition to the requirement of continuity across x = 
x c , Eq. (|4"a|) , we have the same condititions as above apart 
from Eq. (|4ci|) which requires that also the energy current 
is conserved across x = x c . Physically, we assume that 
the superconductor on top of the graphene flake acts as 
an electrode for the high-energy electrons and does not 
noticeably affect the low-energy ones. We also assume 
above and in the following that A is constant: indepen- 
dent of position or heating in the system. 

The degree of inelastic scattering determines the state 
of the system which, in the presence of superconductors, 
can become a quite complicated nonequilibrium state 
when e-e relaxation is incompleted On the other hand, 
when e-e relaxation is complete, two separate alterna- 
tives are possible: The system can either be in the equi- 
librium state fo(e) with bath temperature Tbath or in a 
so-called quasiequilibrium state /o (e) with electron tem- 
perature T e > Tbath depending on whether the region in 
question is heated in one way or another. In our config- 
uration, regions 1 and 2 are separated at x = x c by a 
grounding superconductor effective at e > A. Therefore, 
for the ideal case of complete e-e relaxation, region 2 re- 
mains in equilibrium at temperature Tbath and region 1 
in quasiequilibrium with T e determined by the heating 
voltage V. However, when e-e relaxation is incomplete, 
the presence of the heat link between the regions makes 
it a priori possible that the system in region 2 is in any 
of these states: equilibrium, quasiequilibrium or nonequi- 
librium. Irrespective of the state of the system, we can 
illustrate the heat distribution in the system by defining 
a (local) effective electron temperature 



feeTs = 




where 6(e) is the Heaviside step function. This equates 
the energy in the system to the thermal energy. 

For the quantitative results, we solve the distribution 
function /(e) from the diffusion equation Eq. (flj and 
when V = 0, we use the equilibrium Fermi function 
/(e) = fo(e) with a given bath temperature Tbath- In 
both cases, the supercurrent through graphene is^ 

W) = ^J *is(e,0)[l ~ 2/(e)], (6) 

where we average /(e) over the rr-coordinate. Here, 
C < 1 is a prefactor of the order unity describing the 
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imperfections in the measurement. The spectral super- 
current js{e,4>) °f a diffusive SNS system is defined in 
Ref. [HI and we determine it numerically at the phase 
<t> = 0.67T, which gives a fair approximation for the low- 
temperature critical current I c = maxis(0) ss is(0.67r). 

We note that for graphcne the diffusion constant D 
can depend strongly on the type of disorder and on dop- 
ing. For screened Coulomb impurities ; 14 ' 25 D cx r oc 
Below, D is essentially fit to the experiments, so that we 
do not need to specify the nature of the scatterers. The 
use of a semiclassical diffusion approach has some lim- 
itations for the description of graphene, however. The 
experimentally found finite minimal conductivity at the 
CNP cannot be explained without quantum-mechanical 
effects and a more detailed description of the impurities 
and the nonuniform doping effects they may introduce 
(charge puddles). These are not an issue if the graphcne 
is strongly n or p doped. In our case, for V g > Vcnp the 
doping is of n type everywhere, but as mentioned above, 
for V g < Vcnp p-n junctions are expected to emerge close 
to the contacts^ - — For these reasons, and because of 
the assumed one-dimensionality, our description should 
be viewed only as an effective model, in particular when 
V g < Vcnp- 



B. Inelastic interactions: electron-phonon 

We start by estimating the strength of e-ph contribu- 
tion in the inelastic collision integral. For this, we deem 
it sufficient to pay attention only to acoustic phonons for 
the range of temperatures and voltages relevant to our 
experiments (T,eV/kB ~ 1 K = 0.1 meV). The colli- 
sion integral for acoustic phonons in graphcne has been 
derived in Ref. HI but here we only need the e-ph power 
discussed in Refs. l9l. ITol. l20l andl2ll. We assume the dis- 
tribution /(e) is sufficiently well defined by the effective 
electron temperature T e of Eq. ([5]). Then, in the limit 
(c/vf)\/j,\ ^> ksT e , where c is the sound velocity, the 
c-ph power is 

P c _ ph = S(/i)A(T e 4 - T b 4 ath ), (7) 

with a ^-dependent interaction constant £(/^) and the 
area of the flake A. This power law is applicable up to 
very close to the Dirac point since the condition /i = 
HvFkF = hvpy/ ne r eoSV g Jed kBT e vp/c translates to 

/ kf>T„ \ 2 ed 

5V a » -f-* 8 

\ he J 7re r e 

for the distance from the Dirac point in terms of gate 
voltage SV g = \V g — Vcnp |- For T e ~ 1 K, thickness of the 
gate oxide d = 250 nm, relative permittivity e r = 4 and 
c = Q.Q2v F = 2 x 10 4 m/s, this becomes 6V g > 15 mV. 
In our experimental data the minimal SV g is about 1.5 
V, and so the condition is not violated. This estimate 
neglects the effect of charge puddles close to the CNP, 



but we expect Eq. ((?]) to hold whenever the use of our 
semiclassical approach is justified. 

The total power injected into the system is Pi n = 
V 2 /Rn and if electron-phonon coupling is absent, all of 
this escapes into the leads. If, furthermore, e-e inter- 
actions are neglected, the electron distribution function 
is comprised of discrete steps due to multiple Andrccv 
reflections^ so that highest local effective temperature 
of the system is roughly k B T™ ax « (A + eV)/4M As a 
result, energy escape takes place also at eV <C A. In 
the presence of e-e interactions, T e is reduced from this 
value so that for strong interactions, ksT^ 1 ^ rj eV/A as 
in a diffusive conductor without superconductivity (No- 
tice that, at high voltages, a hot spot is formed in the 
middle of the normal-conducting region 1, so that T™ 3 * 
becomes high while the average T e remains much lower) . 
For T e » T bath 

— (9) 

and we see that the relative importance of e-ph interac- 
tion increases in two cases: first with increasing eV and 
second with eV decreasing below A when e-e interaction 
is absent. We evaluate the both possibilities to obtain a 
range for eV where P c -ph is significant. 

The value for the e-ph interaction constant has 
not been measured but we obtain a theoretical up- 
per estimate S = k%V 2 Tr 2 \fj,\/(60H 5 pv F c 3 ) < 6.9 x 
10~ 3 W/K 4 m 2 for graphene mass density p = 
0.76 mg/m 2 , deformati on p otential constant V = 10 eV 
(see, for example, Ref. |22I ) and chemical potential /i = 
0.22 eV, corresponding to V g = +30 V. This is roughly 
an order of magnitude smaller than e-ph interaction in 
metals j 1 ^ when the reduced dimensionality of graphene 
is taken into account using thickness ~ 1 A. Setting 
A rj 2.8 irni x4.0 /im, A/k B « 1.12 K, and R N w 1 kO, 
yields < 0.013 x (fc B T e /A) 4 /(eV r /A) 2 . We then 

have the result: 

^<"^(£)'(^) 4 - ™ 

assuming ksT e = (A + eV)/4. The ratio is above 1% 
when eV < 0.083A or eV > 12A, while our measure- 
ments are focused on the range eV = 0.5 . . . 5A. We 
emphasize that in this estimate we used the hot-spot 
temperature which is much larger than the average T e 
at high voltages. On the other hand, at low voltages, we 
assumed a total absence of e-e interaction resulting in a 
high T e of the order A. Generally, T e can be expected to 
be even smaller than estimated above and we conclude 
that for our experiment e-ph coupling can be neglected. 

C. Inelastic interactions: electron-electron 

Coulomb interactions in graphene have been discussed, 
for example, in Refs. 0,0, and[l7l However, most of the 
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existing results arc for clean, charge-neutral graphene, 
where the golden-rule collision integrals are furthermore 
plaqued by divergences.— A full theory of e-e interactions 
for diffusive graphene that would be valid at both the 
Dirac point and at finite doping is currently lacking. In 
particular it is not known if a well-defined quasiequi- 
librium state ever exists in diffusive graphene biased 
far from equilibrium, although it is often a convenient 
assumption*^ Since the interband relaxation due to e- 
c collisions is expected to be weak , 6 ' 18 i 19 electrons and 
holes (or electrons in the conduction and valence bands) 
may in any case have to be treated separately^ 

As explained above, our semiclassical approach re- 
stricts our calculation in principle to the strongly-doped 
regime, where only one charge carrier is dominant. In 
this case the system may be expected to behave some- 
what similarly to other disordered two-dimensional con- 
ductors. Thus, a reasonable starting point for an effective 
description is the Altshuler-Aronov theory 2 ^ for diffusive 
normal metals. The collision integral for a well-screened 
diffusive wire in two dimensions i a 26 ' 28 



,(/(<0) = «e 



d(hu) / de'(huj)- 1 



x [I in (uj,e,e')~r ut (uj,e,e')] 



(11) 



with P»(w, e, e') = [1 - /(e)] [1 - f(e')]f(e - foj)/(e' + M 
and 7 out ( W ,e,e') = f(e)f(e')[l-f(e-fuv)}[l-f(e' + hw)}. 
The prefactor is given by ft e _ e = l/(8|/x|r) and while the 
distribution functions arc assumed position-dependent, 
we disregard any such dependence in K e _ e for simplicity. 

We can also estimate the e-e relaxation rate, which 
reada^ 



1 



kBT e . Ti 

In — 

2m P D T„ 



(12) 



in energy-averaged form valid at low-energies ~ fc B T e . 
The result is obtained from the collision integral, where 
an energy cut-off is required to remove a logarithmic di- 
vergence at small u, characteristic to two-dimensional 
systems^ The cut-off energy is Hloq = AksTp/T 2 and 
this gives T x = m 2 e D 3 e 4 N 2 /h 2 k B el, also for Eq. ([12]). 
Here, m e is the electron mass, £o the vacuum permittiv- 
ity and Np the density of states at the Fermi level. For 
graphene, the latter becomes^! 



l/'l 



(13) 



including spin degeneracy. 



D. Final model in dimensionless form 

With the e-e interaction as the sole contributor to the 
inelastic relaxation, we now cast the diffusion equation 
Eq. JI) in a dimensionless form directly applicable for 
numerical solution: 



where x = x/L with L denoting the total length of the 
graphene flake. The dimensionless parameter K e _ e de- 
scribing the strength of e-e interaction in two dimensions 
is derived from K e _ e given above and expressed in terms 
of experimentally relevant parameters. When all energies 
are normalized by A, 
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4 Rq I?Th L ' 



(15) 



where we have used the formula a = e NpD = L/Rpjw 
for conductivity and Eqs. ([2j and (|13p for D and Np- 
Here, Rjy is the normal-state resistance of the graphene 
flake, Rq = h/e 2 the quantum of resistance, i?Th = 
HD/L 2 the Thouless energy, w the width of the flake and 
L its length. Note that K e - e now depends on the dimen- 
sions of the sample and we define it, together with the 
other extensive materials parameters, here for the whole 
flake. Since -Erh oc D and Rn oc D , the parameter 
K e - e is predicted to scale as K e _ e oc D~ 2 oc R 2 N . 

We note again that while K e _ e does not depend explic- 
itly on fi, such dependence is in principle present through 
the diffusion constant D(fi) (and hence Rn and -Erh)- 
Below, we use only the single parameter K e - e to charac- 
terize the e-e interaction and extract it from the experi- 
ments. 



IV. RESULTS 
A. Electron-electron strength 

In order to access the e-e scattering strength K e - e , 
we make two measurements. In the first one, the bath 
temperature is varied. In the second, an injection voltage 
at the left superconductor is used to heat up the graphene 
flake. The critical current and the retrapping current in 
region 2 are then measured as functions of temperature 
and voltage, respectively. In the first measurement, the 
system is in thermal equilibrium, whereas in the second 
one, it can be in a nonequilibrium state. Therefore, we 
use the first measurement to determine i/rh in regions 
1 and 2, and the second measurement to compare the 
experiments with our theoretical model. 

In thermal equilibrium, wc calculate the theoretical 
value for I c simply by using an equilibrium function 
fo at temperature Tb a th in Eq ((6|). Here and below, 
A = 100 /j,eV = 1.1 K so that the critical temperature 
T c ~ 0.6 K. Since js(e, </>) in Eq. ^ is dependent on E^h, 
we may fit the calculated value of I c to the measured one 
at different V g (Figs. [7] and [8]), to obtain the values in 
Table |TJ These values for E^h are subsequently used to 
determine the effective lengths L\ and L2 of regions 1 and 
2, respectively. We assume that the diffusion constant D 
is constant throughout the graphene flake so that 



-K T 

J- »- P — f J- P 



(14) 



Li 
L 2 



'^Th (region 2) 
£rh (region 1) 



(16) 
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TABLE I. Graphene properties determined for different gate 
voltages. The results for £+h and C (Eq. (0) are based on 
critical-current measurements for which we have data for all 
gate voltages apart from V g = —30 V and V g = —3.4 V in 
region 1. 







Region 


1 




Region 2 




Vg 


fj,(meV) 


R N 


.E+h/A 


C 




-E-Th/A 


C 


-30 V 


150 


760 n 


N/A 


N/A 


708 n 


1/3.5 


0.33 


-9.5 V 


42 


1100 Q 


1/6.4 


0.11 


547 n 


1/2.35 


0.42 


-3.4 V 


96 


650 n 


N/A 


N/A 


467 n 


1/2.5 


0.39 


+30 V 


220 


330 n 


1/7.7 


0.18 


337 n 


1/2.7 


0.36 



TABLE II. Parameter values determined for the graphene 
flake. Resistance is calculated as the sum of Rm in region 
1 and 2, whereas Thouless energy is obtained assuming that 
the diffusion constant determined from the J c -measurements 
in region 2 is homogeneous throughout the whole flake. The 
theoretical value for K e - e is calculated from Eq. (|15l) . 



The whole flake 










v a 


Rn 


■Exh/A 


D 




T^theory 

e — e 


-30 V 


1468 Q 


1/24.9 


0.046 m 2 /s 


30.8 


0.51 


-9.5 V 


1647 Q 


1/16.6 


0.069 m 2 /s 


51.9 


0.38 


-3.4 V 


1117 n 


1/17.8 


0.065 m 2 /s 


26.6 


0.28 


+30 V 


667 n 


1/19.5 


0.060 m 2 /s 


7.6 


0.18 




0.1 0.2 0.3 0.1 0.2 0.3 



FIG. 7. (color online). Critical current in region 1 as a func- 
tion of bath temperature. Solid lines are the fits to the theory, 
using £+h and C as the fitting parameters. 
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FIG. 8. (color online). Critical current in region 2 as a func- 
tion of bath temperature. Solid lines are the fits to the theory, 
using Em and C as the fitting parameters. 



This approximate value is used in the simulation and 
assuming L = L\ + L2, we have L\ = 1.75 fivo. and 
L-2 = 1.05 /xm. For consistency, we may also check the 
clastic scattering length I = ILpE^jfiVF resulting from 
the values we obtained for E+h and depending on V g , 
we have I = 90 ... 140 nm. At V g = 30 V, the value 




eV/A 



FIG. 9. (color online). Critical current in region 2 as a func- 
tion of heater voltage. Solid lines are the fits to the theory, 
using K e -e as the fitting parameter at a single point eV < A 
(The exact value depends on the gate voltage) . 



I = 120 nm is relatively close to the experimentally de- 
termined I = 70 nm and the difference can be due to 
inaccuracy in determining the effective length. 

In the second experiment we use the heater voltage 
V at Tbath = 80 mK and, as a result, I c in region 2 is 
dependent on the strength of the e-e interaction K e _ e . 
We proceed by fitting the theoretical and experimental 
results at a single small voltage eV < A with K e _ e as the 
fit parameter. We then use the fitted K e - e to compute 
the whole I c {V)-cmve. The results for different V g are 
given in Fig. [S] The corresponding K e - e are listed in 
Table HH 

From Eq. fTS"]) we expect the e-e strength to scale as 
K e _ e ~ Rjj. The values of K e - e extracted from the ex- 
periments arc shown in Fig. [TO] as a function of Rn from 
Table HT1 assuming additionally that Rn — > 0, K e - e -+ 0. 
The data supports the quadratic behavior. 
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FIG. 10. (color online) . Values for the e-e interaction strength 
K e ~ e determined from the experimental data and shown with 
respect to the normal state resistance of the graphene flake 
Rn- 



B. State of the system: Quasi- or nonequilibrium? 
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FIG. 11. (color online). Electron distribution functions /(e) 
for a representative case (V g = —30 V, eV = 2.5A) at different 
positions of the graphene flake. Here, e-e interactions are 
neglected, i.e., K e - e = 0. Dashed line is the quasiequilibrium 
function /o(e) at the effective temperature T e corresponding 
to the nonequilibrium function /(e). 



Wc determine the state of the system and assess how 
sensitive it is to any changes in K e _ e by looking at the 
electron distribution functions /(e). In our theoretical 
model based on metallic e-e interaction, gate voltage has 
no explicit effect on the results and any gate dependence 
is due only to such dependence in the materials parame- 
ters such as Rn. We therefore use the parameter values 
from the case where 7 ff = -30 V as a representative ex- 
ample. In addition, we set eV = 2.5A so that any pecu- 
larities due to the nonequilibrium state should be visible 
both in region 1 and 2. A true nonequilibrium distribu- 
tion with A e _ e — is given in Fig. llll and we notice that, 
in the absence of relaxation, /(e) remains constant in re- 
gion 2. This happens because the electrostatic potential 
at the middle superconductor is fixed. The numerically 
determined /(e) for AT e _ e = 30.8 (Table HI]) is shown in 
Fig. [12] We see that while the electrons clearly have 
a nonequilibrium distribution near the superconductors, 
where /(e) in any case is strongly affected by the bound- 
ary conditions of Eqs. © and (QJ, in the middle of both 
regions 1 and 2 /(e) is smoothed very close to a thermal 
distribution. 

The question now regarding the state of the system is: 
is /(e) in region 2 effectively a thermal (quasi)cquilibrium 
distribution? If so, the e-e interaction is weak enough to 
let some of the energy injected into region 1 leak into 
region 2, but still so strong that it forces the system 
in region 2 into a thermal state. To answer this, we 
calculate I c function of heater voltage for several 
values of A e _ e in Fig. [13] First, the expected value of 
supercurrent is clearly dependent on the value of A" e _ e in 
this range so the measurement can be used to determine 
K e _ e with satisfactory precision. Second, as there is no 



plateau of constant I c at small voltages, the system is 
only in complete equilibrium state at V = 0. Third, at 
voltages eV > A there is a clear difference between the 
critical current obtained using actual nonequilibrium dis- 
tributions and their quasiequilibrium counterparts, i.e., 
equilibrium functions /o(e) with T e determined from /(e). 
This implies that not only electron heating but also the 
formation of nonequilibrium state affects the observable 
supercurrent at eV > A. For voltages smaller than this, 
a quasiequilibrium description for the system is apt. 



C. Relaxation time 



We finally estimate the e-e relaxation time in our sam- 
ple. From Eq. ([j"2"]h wc find T e _ e = 2 . . . 4 ns at T e = 
Tbath = 80 mK and r e - e = 400 . . . 700 ps at T e = 0.5 K, 
depending on the gate voltage. The relaxation rate is ex- 
pected to be proportional to the electron-electron scat- 
terin strength K e - e . Our measurements therefore sug- 
gest e-e relaxation times which are smaller than the val- 
ues above by a factor K^LJ 'Kf^° ry . Then, at T e = Tbath, 
T e _ e = 30 . . . 70 ps and T e = 0.5 K, r e _ e = 5 . . . 13 ps. We 
are not aware of any comparable results in the low-energy 
regime. The e-e scattering times have been estimated for 
graphene in Ref. [It], but there the graphene is ballistic, 
the energy scale much larger and the obtained scattering 
times, consequently, (even) much smaller. 
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FIG. 12. (color online). Electron distribution functions /(e) 
for a representative case (V a = —30 V, eV = 2.5A) at different 
positions of the graphene flake. Here, the strength of the 
e-e interaction corresponds to the value of Table [TT] which 
we obtained by fitting to the experimental data: K e - e = 
30.8. Dashed line is the quasiequilibrium function /o(e) at the 
effective temperature T e corresponding to the nonequilibrium 
function /(e). The spectral supercurrent js(e) (Eq. (JSJl) is 
shown as a black dotted line with arbitrary units and centered 
so that j s = when |e| < E T h = A/3.5. 
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eV/A 

FIG. 13. (color online). Calculated critical current as a func- 
tion of voltage for different values of the e-e strength K e ~ e . 
The solid lines are the results from the nonequilibrium dis- 
tributions /(e), whereas the dashed lines demonstrate how 
Ic differs if quasiequilibrium functions with effective electron 
temperatures are used instead. Blue crosses correspond to 
the measured values. 



V. DISCUSSION 

We have measured the strength of e-e interaction in 
graphene at four different gate voltages. From Eq. (TT5)) 
we find the expected values for the parameter K e - e de- 
scribing the interaction strength as given in Table HU The 
measured values arc roughly 40-140 times larger than ex- 
pected from the Altshuler-Aronov theory with the largest 
differences closer to the Dirac point. Discrepancies be- 
tween theory and experiments are reported also in metal- 
lic wires (Ag)^ but there the differences are sample- 
specific and only up to a factor of 20, with the measured 
value larger there as well. Even though the strength of 
the e-e interaction is stronger than expected, the system 
is still not thermalized and the incomplete e-e relaxation 
can be seen at heater voltages V well below the supercon- 
ducting gap by measuring the critical current. The result 
at low voltages is seen as an energy leak from the heater 
junction and increasing temperature in the thermometer 
region. For eV > A, also the electrons in the thermome- 
ter are driven into a nonequilibrium state where a ther- 
mal description with an effective temperature T e is no 
longer enough. In addition to finding the magnitude of 
the e-e scattering strength, wc find that the scattering 
strength exhibits a significant gate dependence, presum- 
ably due to changes in charge density as the gate voltage 
is varied. 

Finally, we note that we have also estimated the inter- 
action strength between electrons and acoustic phonos in 
our setup with an aim to measure it. However, since the 
e-e interaction is relatively strong, we expect that T e is 
at most of the order of the heater voltage except when 
eV <C A. Consequently, the expected electron-phonon 
power at eV < A becomes even lower than predicted by 
our estimate in Sec. IIII Bl At the other end of the scale, 
eV > A, the effective temperature needs to be increased 
yet further for a high ratio P c - P h/-Pin ~ T^/V 2 . This re- 
sults in notable heating also in the thermometer region, 
making a critical current measurement such as the one 
used here very difficult, unless a strong thermal isolation 
is established between the heater and the thermometer 
regions. 
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